Angular dependence of Josephson currents in unconventional superconducting junctions 
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Josephson effect in junctions between unconventional superconductors is studied theoretically within the 
model describing the effects of interface roughness. The particularly important issue of applicability of the 
frequently used Sigrist-Rice formula for Josephson current in d-wave superconductor / insulator / d-wave super- 
conductor junctions is addressed. We show that although the SR formula is not applicable in the ballistic case, 
it works well for rough interfaces when the diffusive normal metal regions exist between the d-wave supercon- 
ductor and the insulator. It is shown that the SR approach only takes into account the component of the d-wave 
pair potential symmetric with respect to an inversion around the plane perpendicular to the interface. Similar 
formula can be derived for general unconventional superconductors with arbitrary angular momentum /. 
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A number of phase- sensitive experiments have convinc- 
ingly demonstrated the realization of d-wave pairing state in 
high Tc cuprates QESIil]- Because of such unconventional 
symmetry, the study of Josephson effect in high Tq supercon- 
ducting (HTS) junctions attracted a lot of interest. A while 
ago simple formula for the Josephson current of d-wave su- 
perconductor / insulator / d-wave superconductor (DID) junc- 
tions was proposed by Sigrist and Rice (SR)|5]. According 
to the SR formula, the Josephson current is proportional to 
cos 2 a cos 2(3, where the a(/3) denotes the angle between the 
normal to the interface and the crystal axis of the left(right) 
(i-wave superconductor |5|]. 

Although the SR formula can explain experiments with 
the so-called 7r -junctions |ll I2L this formula does not take 
into account the effect of mid-gap Andreev resonant states 
(MARS) formed at junction interfaces 0, 01- Actually, as 
shown in Ref. 8 , SR formula does not work in ballistic d- 
wave junctions for a ^ and (3 ^ where MARS influ- 
ence severely the charge transport at low temperatures. It was 
shown both theoretically |9, 10] and experimentally II llll2ll 
that MARS induce a nonmonotonic temperature dependence 
of the maximum Josephson current in DID junctions. On the 
other hand, SR formula has been extensively used to analyze 
experiments with various types of HTS Josephson junctions 
L^UilLLj] Experiments with HTS junctions are of high im- 
portance for basic understanding of high-Tc superconductiv- 
ity since they may provide an information on poss ible sub- 
dominant admixtures to the d-wave symmetry |16l [nl El- 
Therefore it is of fundamental interest to understand the phys- 
ical mechanisms which determine the angular dependence of 
Josephson current in HTS junctions. For this reason, the deter- 
mination of the conditions of applicability of the SR formula 
is an important issue which is addressed in the present paper. 

In the following, we study the Josephson current in 
D/DN/I/DN/D junctions, where DN denotes diffusive normal 
metal and could be formed between the insulator and d-wave 
superconductors. The calculations are based on the quasiclas- 
sical Green's function method applicable to unconventional 



superconductor junctions We find that the resulting 

Josephson current in D/DN/I/DN/D junctions is well fitted by 
the SR formula. Near the transition temperature, it is proven 
analytically that Josephson current follows the SR formula. 
We also confirm that this formula does not hold in the ballistic 
junctions. It is clarified that in the SR formula, the component 
of the pair potential which is antisymmetric by the inversion 
operation around the plane perpendicular to the interface is 
neglected. We also study p-wave superconductor / diffusive 
normal metal /insulator/ diffusive normal metal/ p-wave su- 
perconductor (P/DN/I/DN/P) junctions. The resulting Joseph- 
son current is also well fitted by cos a cos (3, where a {(3) de- 
notes the angle between the crystal axis of left (right) p-wave 
superconductor and the normal to the interface. This is a cor- 
responding version of the SR formula in the p-wave supercon- 
ductor junctions. Furthermore, it is possible to extend the the- 
ory for an unconventional superconductor (US) with arbitrary 
angular momentum I. For US/DN/I/DN/US junctions, the ex- 
pected Josephson current is proportional to cos la cos 1(3. The 
obtained results may serve as a guide for the analysis of the 
experiments in unconventional superconductor junctions. 

Before we proceed with a formal discussion, let us pro- 
vide qualitative arguments on the physical meaning of the 
SR formula and explain why it holds in the diffusive junc- 
tions. First, we consider d-wave superconductor junctions. 
The pair potentials of left and right d-wave superconduc- 
tors are, respectively, expressed by A^ = A[/sx(0) + 
/asl(0)]bxpH*), and A R = A[f SR ^) + Iasr^)}, 
with fsL,(<t>) ~ cos 20 cos 2a , /asl(0) = sin 20 sin 2a, 
Ssr{§) = cos 20 cos 2/3, fASR{4>) = sin 20 sin 2(3 , where 
is the injection angle measured from the interface normal, 
A denotes the maximum value of the pair potential and \I/ 
is the phase difference across the junction. The terms pro- 
portional to cos 20, i.e., /5l(0) and /s j r(0), correspond 
to the d x 2_ ?/ 2-wave pair potential and the the terms propor- 
tional to sin (20), i.e., /asl(0) and /asrW, correspond to 
the d X y-wave pair potential, respectively. Here, /sx(0) = 

/5L(-0), fsR{4>) = fsR{-4>), fASL(4>) = ~/a5l(-0), 
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and fASR(<t>) = —fASR(—^) are satisfied. In the actual cal- 
culation of Josephson current in D/DN/I/DN/D junctions, we 
have to take an average over the various cj). Due to the impu- 
rity scattering in DN, the average is taken for the left and right 
D/DN interface independently. Then we can drop Jasl^) 
and /asrW and arrive at the SR formula, where only the 
terms fsL,{4>) and fsR{4>) remain which do not change sign 
by exchanging for —cj). This fact is in accordance with the 
recent result that the proximity effect is absent in the case 
of (i^-wave pair potential [19]. Consequently, the resulting 
Josephson current is proportional to cos 2a cos 2/3. 

Similar arguments apply to p-wave junctions. In this case, 
fsL(4>) = cos cos a, fASL(4>) = sin sin a, fsR^) = 
cos <j) cos (3, and /asrW = sin^sin/3 are satisfied. The 
terms proportional to cos0, i.e., fsL,(<j>) and fsR((/>), corre- 
spond to the p^-wave pair potential and the the terms pro- 
portional to shi(/>, i.e., /a5l(0) and /asrW, correspond 
to the p^-wave pair potential. In the actual calculation for 
P/DN/I/DN/P junctions, functions /aslW and /asrW 
vanish after averaging over angle <p. This is consistent with 
our previous results that the pair potential with p y -wave sym- 
metry does not contribute to the proximity effect 1EHE2I1 . 

Next we formulate the junction model and basic equations 
starting from the <i-wave case. We consider ballistic DID 
and D/DN/I/DN/D junctions. The DN has a resistance R d 
and a length L much larger than the mean free path. The 
DN/D interfaces located at x = ±L have the resistance R' h , 
while the DN/I interface at x = has the resistance R^. We 
model infinitely narrow insulating barriers by the delta func- 
tion U(x) = H / 5(x^L)^HS(x)^H / S(x-L). The resulting 
transparencies of the junctions T m and T' m are given by T m — 
4 cos 2 0/(4 cos 2 0+Z 2 ) and = 4 cos 2 0/(4 cos 2 (j)+Z' 2 ), 
where Z = 2H/vp and Z' = 2H' jvp are dimension- 
less constants and vp is Fermi velocity. Below we assume 
Z ^> 1. The schematic illustration of the models is shown 
in Fig. [I] Here, a and (3 denote the angles between the nor- 
mal to the interface and the crystal axis of the left and right 
(i-wave (or p-wave) superconductors, respectively. The lobe 
direction of the pair potential and the direction of the crystal 
axis are chosen to be the same. The pair potential along the 
quasiparticle trajectory with the injection angle </> is given by 
A L = Acos[2(0-a)] exp(-i*) and = A cos [2 (</>-/?)] 
for the left and the right superconductor, respectively. For bal- 
listic junctions, we use a similar model without DN and cal- 
culate the Josephson current following Ref. [9]. 

We parameterize the quasiclassical Green's functions G and 
F with a function lEiEl: 
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where w is the Matsubara frequency. Then the Usadel equa- 
tion reads 1 25] 
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FIG. 1: (color online) Schematic illustration of the models of (a) 
D/DN/I/DN/D and (b) P/DN/I/DN/P junctions. 



with the coherence length £ = ^/D/2ttTc, the diffusion 
constant D and the transition temperature Tq- We solve the 
Usadel equation with the boundary conditions in Ref. lEoll at 
x = ±L and those in Ref J26l1 at x = 0. 
The Josephson current is given by 



eIR 
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where T is temperature and R = 2R d + R^ + 2R' h is the 
normal state resistance of the junction. In the following we 
focus on the IcR value where Iq denotes the magnitude 
of the maximum Josephson current. We fix parameters as 
Z' = 0, R d /R b = 0.01, R d /R' b = 10 and E Th /A = 1 
for D/DN/I/DN/D junctions and Z = 10 for DID junctions. 
Ao denotes the value of A at zero temperature. The choice of 
the small magnitude of Z 1 and R' h and large Thouless energy 
is justified by the fact that thin DN is naturally formed due to 
the degradation of superconductivity near the interface. 

The a dependence of IcR for (i-wave superconductor junc- 
tions is plotted in Fig. |2] In Figs. |2fa) and |2fb), IcR of 
ballistic junctions is plotted for low (T/Tc = 0.2) and high 
temperature (T/Tc = 0.9), respectively. With the increase 
of the magnitude of f3, the dependence of IcR on a trans- 
forms from cos 2a to sin 2a. These a dependences can not be 
expressed by the SR formula, where IcR is proportional to 
cos 2a for fixed (3. On the other hand in D/DN/I/DN/D junc- 
tions, IcR has a simple form, cos 2a, independent of (3 at low 
and high temperatures as shown in Figs.|2tc) and (d), respec- 
tively. The magnitudes of IcR in D/DN/I/DN/D junctions are 
at least two orders smaller than those in DID junctions. By 
taking account of the (3 dependence, IqR is almost propor- 
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tional to cos 2a cos 2(3. It should be remarked that this fitting 
is possible for small magnitude of Z' and R' h /Rd where the 
MARS formed at the D/DN interface do not influence seri- 
ously the charge transport. 

The corresponding results of IcR for p- wave superconduc- 
tor junctions are plotted in Fig. [3] For P/DN/I/DN/P junc- 
tions, IcR can be fitted by cos a cos (3. Similar to the case of 
(i-wave junctions, this fitting is possible for small magnitude 
of Z' and R' h /Rd. For ballistic junctions, as shown in Figs. 
|3a) and[3b), this fitting does not work any more. 
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tions and works in the diffusive junctions. Although we focus 
on d- wave junctions, similar discussion is possible for p- wave 
junctions. Near Tq (A <C u), we can get the formula for the 
ballistic DID junctions 1 9] : 



eIR A 2 sin ^ 



(4) 



ttT c STT c 
-< cos 2 2(j) > cos 2a cos 2/3+ < sin 2 20 > sin 2a sin 2/3 



Here, the average over the various angles of injected particles 
at the interfaces is defined as 



< B{<t>) >-- 



/:f /2 #T(0) cos 05(0) 



/^/ 2 #T(0) COS0 



(5) 



with T(0) = T m . It is easy to check that (cos 2 20) and 
(sin 2 20} are of the same order for all Z. Therefore, SR for- 
mula cannot be applicable to nonzero values of a and (3. Also 
we can roughly estimate the Josephson current: 
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which is consistent with the result in Fig. E|(b). In the case 
of PIP junctions, we can obtain the corresponding equation 
by replacing 2a, 2(3 and 20 with a, (3, and in the above 
equations. 



FIG. 2: (color online) Maximum Josephson current for d- wave junc- 
tions, (a) and (b) DID junctions, (c) and (d) D/DN/I/DN/D junctions. 
Solid lines in (c,d) are proportional to cos 2a cos 2/3. 
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FIG. 3: (color online) Maximum Josephson current for p- wave junc- 
tions.^) and (b) PIP junctions, (c) and (d) P/DN/I/DN/P junctions 
with solid lines which are proportional to cos a cos (3. 

In the following, we will present analytical result demon- 
strating why the SR formula does not work in ballistic junc- 
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FIG. 4: (color online) Schematic illustration of the inversion sym- 
metry of the pair potential around the plane normal to the interface. 
(a)c/ cc 2_ y 2-wave, (b)p x -wave, (c)d xy -wave, and (d)p y -wave. The 
dependences are given by cos 20, sin 20, cos 0, and sin 0. 

Next we consider the D/DN/I/DN/D junctions. Near Tc, 
we can linearize the Usadel equation as follows, 



^2 ® 2 * 00 ^ 







(7) 



where j(=l, 2) denotes the left or right DN. Similarly the 
boundary conditions at £ = —L,x — and x = L are re- 
duced to 
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-^L \ x=L= _A ( _$ 2u + Iq cos 2(3) \ X=L (10) 

with 7 = A (cos 20). 

Solving the above equations, we find the expression for the 
Josephson current of the form 

eIR Rrr' 2 T ^ jL < cos 2cj) > 2 A 2 cos 2a cos 2(3 sin \£ 



lu 2 F 1 F 2 



(ID 



Fi = 7L sinh7L + r COSI17L, 



F 2 = [(2rr + 7 2 £ 2 ) sinh 7L + (2r + r'^L cosh jL] 



with r = ^, r' and 7 = yJ4*. Thus the SR formula 

is proven to be valid near T c . In the case of P/DN/I/DN//P 
junctions, cos 20, cos 2a and cos 2/3 have to be replaced with 
cos 0, cos a and cos (3, respectively to obtain the correspond- 
ing formula. This result is consistent with the previous study 
of DID junctions with rough interface [ 27] where the SR for- 
mula is applicable as well. 

In order to understand the above results qualitatively, let's 
discuss the symmetry of the pair potential by the inversion 
operation around the plane perpendicular to the interface. As 
shown in Fig. 4, ^2 ^2 -wave and p^-wave are symmetric 
while d^-wave and p^-wave are antisymmetric by this oper- 
ation. Only the symmetric pair wave function is taken into 
account in the SR formula. Applying this idea to an arbi- 
trary unconventional superconductor with angular momentum 
I, one can argue that the Josephson current is proportional to 
cos la cos 1(3. It is straightforward to get this result just by 
replacing cos 20, cos 2a and cos 2(3 with cos/0, cos la and 
cos 1(3 in Eq. (11), respectively. 

In summary, we have studied the validity and the physi- 
cal meaning of the Sigrist-Rice formula in d-wave supercon- 
ductor junctions. According to the SR formula, the ampli- 
tude of the maximum Josephson current is proportional to 
cos 2a cos 2(3. Although this formula is not applicable to the 
ballistic junctions, it works well for D/DN/I/DN/D junctions 
where the DN regions are located between the d-wave su- 
perconductor and the insulator. We have also shown that in 
P/DN/I/DN/P junctions, the Josephson current is proportional 
to cos a cos (3. The obtained results may help to obtain infor- 
mation about pairing symmetry in experiments with uncon- 
ventional superconducting junctions. 
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